Intensity-dependent dispersion under conditions of electromagnetically induced 
transparency in coherently prepared multi-state atoms. 
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Interest in lossless nonlinearities has focussed on the the dispersive properties of A systems under 
conditions of electromagnetically induced transparency (EIT). We generalize the lambda system 
by introducing further degenerate states to realize a 'Chain A' atom where multiple coupling of the 
probe field significantly enhances the intensity dependent dispersion without compromising the EIT 
condition. 

PACS numbers: 42.50.Gy, 42.50.Hz, 32.80.-t, 42.65.-k 



There has been much interest given lately to the en- 
hancement of optical nonlinearities in Electromagneti- 
cally Induced Transparency (EIT). Most of the work 
has focussed on the three-state system in the A config- 
uration, which has provided some dramatic examples of 
nonlinear optical effects. Examples include ultra-slow 
, stopped j^J and superluminal (^] group velocities, co- 
herent sideband generation Q etc. All these nonlinear 
processes depend on the creation of coherent superpo- 
sitions of the ground states with accompanying loss of 
absorption, and such mechanisms were described in [^. 
Thorough reviews of EIT and its properties can be found 
in and 0. 

Recent investigations of nonlinear optics at the few or 
single photon levels have identified four state systems 
where the probe field simultaneously couples two tran- 
sitions in the N configuration. Examples of applications 
for such work includephoton blockade Q and two photon 
absorptive switches Q . The classical precursors to such 
experiments have also been performed 10] llj. Other 
experiments on the N scheme have been performed by 
Entin et al. [l^ . In order to realize larger nonlinear 
effects, Zubairy et al. ^3 suggested an extension where 
the more usual N configuration was extended to a system 
with an arbitrary (even) number of states where all the 
states are resonantly coupled except on the final tran- 
sition where detuning is present. This scheme shows 
enhanced nonlinearities of not only but also higher or- 
der susceptibilities. One problem with this scheme and 
the standard N scheme is to do with the need to balance 
the required nonlinearity and decoherence in the system. 
To enhance the nonlinearity present in the systems it 
is important for the detuning of the final probe field to 



be minimized, however decreasing the detuning increases 
the amount of the final excited state which is mixed into 
the coherent superposition state, resulting in an increase 
in decoherence and optical losses. This problem is to 
some extent circumvented by the absorptive switch of 
Harris and Yamomoto ^ by exploiting such losses, and 
in photon blockade by using the cavity nonlinearity to 
prevent absorption of the final photon. Still the increase 
in decoherence proves to be a difficulty in experimen- 
tal precursors to these processes and causes problems in 
travelling wave configurations. 

An alternative multistate configuration for investigat- 
ing EIT enhanced nonlinearities is the tripod configura- 
tion, studied recently by Paspalakis and Knight [Ijl and 
earlier considered by Morris and Shore 15]. This sys- 
tem has many of the advantages of the N system, but 
by maintaining superposition states of the three ground 
states also avoids the excess decoherence of the N sys- 
tem. Morris and Shore 15] also mentioned a multi-leg 
extension of the tripod scheme. 

Here we present an alternative extension to the stan- 
dard A configuration which we term as the Chain A con- 
figuration, depicted in Fig. 1. We start with a A atom 
[Fig. 1(a)] with ground states l^i) and I52) and excited 
state ]ei). The gi — ei (32 — ei)transition is excited by a 
probe (coupling) field of frequency LUp (wc) and detuning 



UJr, 



^eigi (Ac 



j). The probe (cou- 



pling) Rabi frequency is P (C). We assume for conve- 
nience that the transition frequencies Weigi and i^exg2 are 
equal and that appropriate selection rules ensure that the 
probe and coupling fields only interact with their des- 
ignated states (a concrete example of how this can be 
achieved is described below). 
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FIG. 1: Energy level configurations for Chain A atoms, (a) 
is the usual A system, (b) is the 5 state Chain A or M system, 
(c) shows the generalization to higher number of states. 



The five-state Chain A is an M system and is illus- 
trated in Fig. 1(b). Here the ground states are labelled 
\gi), |(72), I53), and the excited states |ei), |e2). The 
probe field excites the gi — ei and (72 — 62 transitions 
simultaneously with Rabi frequencies Pi and P2 respec- 
tively (different labels are applied to take account of the 
different coupling strengths of the transitions), whilst the 
coupling field excites the g2 — ei and g^ — 62 transitions 
with Rabi frequencies Ci and C2. We note that a com- 
plementary work by Matsko et al. 16], which looks at 
Faraday rotation and Kerr nonlinearities in the A, N and 
M schemes, has been performed which confirms some of 
our predictions about nonlinearities in these systems. An 
early study of the M scheme in the context of degener- 
ate two level systems was also performed by Morris and 
Shore [11 



The Chain Lambda atom with 2n — 1 states is shown 
in Fig. 1(c). The ground states are denoted \gi), 
\gn), and the excited states |ei), |e„_i). The probe 
field excites the gj — ej transition with Rabi frequency 
Pj and the coupling field excites the gj+i — Cj transition 
with Rabi frequency Cj. 



In order to gain insight into the problem we consider 
first the Hamiltonian for the Chain A atom with 2n — 1 
states. Using the rotating wave approximation this can 



be written 

n 

h 



^(j-l)(Ap-A,)|5,)(.g,| 

n-l 

J = l 
ra-1 
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We ignore the decay from the excited states in our ana- 
lytical analysis in order to gain simple expressions for the 
dressed states of the field-atoms system and thus gain a 
clearer understanding of the problem. Furthermore we 
shall concentrate our analysis on the optical nonlineari- 
ties which are present in the vicinity of the dark state, 
which are relatively insensitive to decay. This is evident 
by direct comparisons between numerical solutions of the 
complete master equation with decay, and our decay-free 
analytic expressions. The Hamiltonian can be conve- 
niently expressed as a tridiagonal matrix with state or- 
dering |ei), 1^2), \gn), 
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where we have introduced Acp — Ap — Ac. 

Following the approach taken by Kuang et al. [l?| and 
Zubairy et al. [l3j, we first calculate the eigenvectors of 
the Hamiltonian. These can be written as 



I'Di) = ai,3i \gi) + ai,ei |ei) + ai^g^ \g2) H hi 



where i varies from 1 to (2n — 1). It is clearly not pos- 
sible to give general solutions for the a^'s for atoms with 
more than 3 states, although one may simply derive nu- 
merical results. However, if we invoke the adiabatic 
hypothesis [l7j | and assume that the probe detuning is 
small (Ap <C Pi,Ci), the probe field is turned on slowly, 
and the coupling field resonant, then we may assume that 
the system evolves solely into the dressed state with en- 
ergy closest to 0. For convenience we denote this state 

Pq"'^ where n is the number of states in the Chain A 

atom. Using MAPLE [ll|, and the simplification that 
Pi = P, Ci = C, i — 1, 2 ... 71 we have derived expressions 



for 



in) 



These are presented in table as unnormal- 
ized quantities and where fi^ = + P^. Note that the 
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TABLE I: Unnormalized coefficients of \T>o) for 3, 5 and 7 
state Chain A atoms. 
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FIG. 2: A possible way to realize the M system using an 
F = 2to F' — 2 transition. The coupling field is cr^ polarized 
and turned on before the polarized probe. In this way 
the desired M system (bold) is favored over the W system 
(dashed). 



results for 



(3)\ 



^ are equivalent to those which appear 

in Ref . 0| and were also used in Ref . . 

Short of directly creating an artificial atomic structure 
in a quantum well type material, or in an optical lattice, 
it is important to investigate whether the required Chain 
A structure is naturally present in any materials. A sim- 
ple approximate realization to the Chain A configuration 
is obtained by exciting an = n to i^' = n transition in 
an atomic vapor where, as before, 2n — 1 is the number of 
states in the Chain A. We illustrate this for the five-state 
Chain A atom in Fig. 2, although the scheme generalizes 
in an obvious manner. We assume that the coupling 
field is a~ polarized (and therefore excites transitions 
from rriz — n to niz — n ~ 1) and the probe cr"*" polar- 
ized (exciting transitions from = n to = n + 1). 
Notice that without any preparation there are two sys- 
tems here, first the required M system (bold lines in Fig. 
2), and second an undesired W system (dashed lines). 
In order to select the M over the W we first apply the 
coupling field, this has the effect of optically pumping 
the population into the ruz = —2 state. Next the probe 
beam is turned on sufficiently slowly to ensure that the 
system evolves adiabatically to the desired dark state. 
We also note that the W system does not have any dark 
states. So even without the adiabatic state preparation, 
any population in the W system will eventually be opti- 
cally pumped into the desired dark state of the M system. 
We also note that the presence of M systems have been 
identified in conjunction with A systems at least twice 
before in Refs. jl^ It is important to notice that 

our method of realizing Chain A systems and performing 
experiments with them, is not significantly more complex 
than standard experiments on simple A systems, all that 
is necessary to achieve the enhanced nonlinearities is the 
appropriate choice of transition. 

To study intensity-dependent dispersion it is necessary 
to extract the susceptibility at the probe frequency as a 
function of small probe detuning 



X 



i=i ' j=i 
where k = 2TTj\f/eofi, J\f being the atomic density and 



fig e the dipole moment of the gj — Cj transition. The * 
denotes complex conjugation. In order to calculate real- 
istic values of the dispersion which would be attainable 
in standard experiments with alkali atoms (for example 
in a magneto-optical trap or vapor cell), we have taken 
7V= 3xl0iW^ fig^e, = 2xlO-29Cmandr = 5.6MHz. 

In order to derive simple results for the nonlinear dis- 
persions, we shall assume the coupling constants for all 
transitions to be equal, i.e. Hg^ei = M ^^^d so the probe 
(coupling) field Rabi frequencies are the same for all tran- 
sitions Pi ^ P [Ci = C). Using MAPLE we can then 
derive expressions for the intensity dependent dispersion, 
R, for Chain A atoms of varying number of states: 



C72 



1 d 

PdAr, 



E 



We note that our results for i?^^^ are compatible with 
the intensity dependent group velocities derived in Ref. 
[l7j . The analytically determined dispersions are plotted 
in Fig. 3 as a function of P/F for C = 0.25F. In Fig. 
4 we present a comparison between the analytical results 
for the intensity dependent dispersion and that obtained 
by solving the density matrix equation with decay given 
in the Appendix. The good agreement shows that our 
analytic approach is justified. 

If one simply performs a Taylor series expansion on our 
results for the P's (which are non-perturbative), it is easy 
to show that the Chain A systems exhibit nonlinearities 
to all orders in P/C. The linear dispersion has previously 
been identified as being iinportant in KIT systems (see 
for example Harris et al. [13)) and it is clear from Fig. 
3 that the linear dispersion is identical for all Chain A 
atoms. From this we may conclude that in the limit of 
weak probe fields, the probe field cannot couple the levels 
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FIG. 3: Graphs showing dispersion (times hnewidth 
squared), RV'^ as a function of P/T with C/F = 0.25 for 
Chain A systems of 3,5,7 and 9 states. 



FIG. 5: Graphs showing log^p of group velocity, Vg as a 
function of P/T with C/F = 0.25 for Chain A systems of 
3,5,7 and 9 states. 



20 




Q\ , , , , , , , , , 1 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



P/r 

FIG. 4: Graph showing the comparison between intensity de- 
pendent dispersion calculated using the full master equation 
(solid line) and the analytical approach using dressed states 
(dashed line) in a five-state system. 

in any fashion other than the simple A scheme. However 
as the probe intensity is increased, higher-order process 
begin to turn on. In order to rigorously determine the or- 
der of the nonlincaritics present, one should construct an 
effective Hamiltonian, following methods presented by, 
for example, Zubairy et al. jl3l | or Klimov et al. [20j. 
This has not yet been performed for the Chain A system 
and it is hoped that such investigations will shed more 
light on the nonlinear optical properties of these systems. 

An interesting feature to note in the dispersion cal- 
culations is the Rabi frequency ratio which provides the 
maximum dispersion. In this simple case, the position of 
the maximum depends only upon the ratio P/ C . If /? is 



the value of this ratio at the maximum then for P/C < (3 
the dispersion will be monotonically increasing with in- 
creasing P, and for P/C > f3 it monotonically decreases. 
It is clear that reciprocal results will be obtained for the 
corresponding group velocities. For 3, 5, 7 and 9 state 
atoms, the values of (3 are 0, 0.476, 0.698 and 0.804 (to 3 
significant figures) respectively. 

One material property dependent on the dispersion is 
the group velocity, which is 

c 

where rj — ^1 -1- x is the complex refractive index. It is 
important to realize that in systems with nonlinear dis- 
persions as large as those for Chain A systems, the group 
velocity may be a poor parameter. This is because for 
realistic propagation through an optically thick medium, 
the intensity dependence of the medium will alter the 
shape of a simple gaussian pulse. There are however, 
other experiments sensitive to the group velocity which 
may be considered, for example bichromatic excitation of 
the probe beam to generate a beat note or the use of 
a frequency modulated (rather than the more usual am- 
plitude modulated) probe signal. In Fig. 5 we present 
intensity dependent group velocities for 3, 5, 7 and 9 state 
Chain A atoms corresponding to the dispersion calcula- 
tions presented in Fig. 3. 

The analytical results given above cannot provide a 
complete description of the physically realizable prob- 
lem, because of the different Clebsch-Gordan coupling 
between the states involved in the transitions. If we use 
the scheme suggested in Fig. 2 for the couplings and de- 
fine our coupling strengths relative to the coupling in the 
first A system (i.e. 5i — ei — ^2) then we may write down 
the expressions for the other couplings in terms of these 
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TABLE II: Relevant Rabi frequency ratios for 5 and 7 state 
Chain A systems 



stand propagation effects and especially the effect of such 
high nonlinearities on pulse shape. Such analysis goes 
beyond our simple picture and will be the focus of future 
work. 
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FIG. 6: Graphs of intensity dependent dispersions as a 
function of Pi, calculated using realistic ratios of coupling 
strengths, i.e. with the relative Pi and Ci's calculated using 
the Clebsch-Gordan ratios in Table II. 



quantities [2l|. These are summarized in table Hll for the 
5 and 7 state Chain A atoms, and it is easy to generalize 
to higher orders. 

Results of numerical calculations of intensity- 
dependent dispersions for 5- and 7-state Chain A 
systems are presented in Fig. 6. Comparing Figs. 3 
and 6, shows that despite the differences in the values 
obtained for the dispersions, there is only minimal 
change to the overall shape of the intensity dependent 
curves. 

We have shown that there exist interesting nonlinear 
properties for Chain A atoms, and in particular we have 
focussed on the intensity-dependent dispersion as a mea- 
sure for these nonlinear optical properties. The nonlin- 
earity of these systems increases as the number of atomic 
states increases, whilst the FIT transparency is main- 
tained. It therefore appears likely that such multistate 
systems will be useful in the search for new quantum 
non-linear optical materials. Our analysis has been con- 
fined to the optically thin regime. Clearly, in a full study, 
which would include intensity-dependent group velocities 
in such highly nonlinear media, it is important to under- 
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II. APPENDIX 



In this appendix we present the density matrix equa- 
tions of motion for the 5-state Chain A atom (M scheme). 
The equations to be solved are: 



P = --[n,p] + c, 

1 Pgigi ^ Pg2g2 + Pgsga Peiei + Pe 



For reasons of space we split the Hamiltonian superoper- 
ator into smaller sub-blocks, thus 
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where H2j,2j+i = H2j-i,2j 
Tiji- I (5) is the 5x5 identity matrix and 
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FIG. 7: Three-dimensional plot showing — 9[xl (proportional 
to the probe absorption) as a function of Ap/F and P/V for 
Pi = P2 = P, Ci = C2 = C = F/4, F = 5.6MHz for the 
five-state Chain A system. 



for Pi = P2 = P, Ci = C2 = C = r/4, r = 5.6MHz 
and other parameters as above, for the five-state Chain 
A system is presented in Fig. 7. 



where Acp = Ap — Ac. The loss operator is 
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where the off-diagonal blocks {i ^ j) have all ele- 



ments zero except the (ij)*'^ which is r/2. The >Cii's 
are all diagonal matrices, with non-zero elements 

£11 = (0,-r/2,-r2,-r3,-r4), 

£22 = (-r/2,-r,-r/2,-r,-r3), 

Hzz = (Fa, -r/2, 0, -r/2, -Fa), 

A4 = (-r3,-r,-r/2,-r,-r/2), 
^5 = (-r4,-r3,-r2,-r/2,o), 

where r„ is the n photon dephasing and T = Vf.^ = V,.^ 
is the total decay rate from either excited state. 

A three-dimensional plot showing — 3(x) a-s a function 
of probe detuning, Ap/r and probe Rabi frequency, PjV 
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